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THE NEAR-EARTH METEOROID ENVIRONMENT 


By Robert J. Naumann 


SUMMARY 


A detailed study of the meteoroid environment has been completed uti- 
lizing observational data taken from Explorer XVI, Explorer XXIII, Pegasus I, 

II, and III, as well as from radar and photographic techniques. These obser- 
vations cover the entire spectrum of masses important to the meteoroid hazard 
problem. The consistency of these measurements is examined and a mass 
distribution function is found which gives a best fit to the measurements. Pre- 
dictions of penetration frequency are made independently of any mass determi- 
nation which avoids this source of uncertainty. 

The method employed is a differential technique that makes use of two 
or more similar observations with different thresholds to obtain values for the 
logarithmic derivative d log N/d log m, where N is the number of encounters 
per unit area-time solid angle of meteoroids with mass equal to or greater than 
mass m. Since it is only necessary to determine a ratio of masses characteristic 
to different thresholds of a given measurement technique, any multiplicative 
factors relating mass to the observed quantity are cancelled. A differential 
equation derived from the measured log derivatives gives a relation of log m 
and log N to within a single constant. This constant is determined to give a 
best fit of the resulting values estimated from each measurement. A similar 
differential equation is solved to give a relation between log N and log of pene- 
tration to within a single constant which is determined from an appropriate pene- 
tration measurement. 

The mass distribution function was obtained from the cumulative mass 
flux and is found to agree remarkably well with the interplanetary dust distri- 
bution derived by Van de Hulst from interpretation of photometric measurements 
of F-corona and zodiacal light. It was found that interplanetary measurements 
of dust by Alexander on Mariner IV were in order of magnitude agreement. These 
results suggest that the environment is not substantially different in interplanetary 
space. Also, the fact that the slopes determined from satellite penetration data 



appear to rapidly approach zero as N increases, suggests an upper limit on 
number of meteoroid encounters which appear to be on the order of 3/m 2 day 
regardless of size. This low value for the maximum encounter frequency de- 
mands that a meaningful meteoroid experiment must expose an area on the order 
of a m 2 to collect a good statistical sample in a reasonable time. 

The frequency of penetrations predicted from the analysis is much 
smaller for very thin material (< 0. 0025-cm Al) than any previous estimate. 
However, the number of penetrations predicted for materials > 0. 04-cm Al is 
more nearly approximated by the more pessimistic previous estimates. There 
is little doubt that the meteoroid hazard is indeed significant for vehicles with 
large exposures (area-time product) and that considerable emphasis must be 
placed on meteoroid protection in the design of such vehicles. 


INTRODUCTION AND APPROACH 


A detailed study of the meteoroid environment has been completed uti- 
lizing observational data taken from Explorer XVI, Explorer XXIH, Pegasus I, 

H, and III, as well as from radar, and photographic techniques. These obser- 
vations cover the entire spectrum of masses important to the meteoroid hazard 
problem. The consistency of these measurements is examined and a mass dis- 
tribution function is found which gives a best fit to the measurements. Predic- 
tions of penetration frequency are made independently of any mass determination 
which avoids this source of uncertainty. 

The usual approach for obtaining a meteoroid hazard model first attempts 
to obtain a mass flux relation by interpreting various observations in terms of 
mass. Since there is no direct means for measuring mass of a meteoroid, such 
information must be inferred from observables such as penetration in known 
material, electron trail density, or emission of light. Such observables are also 
strongly dependent on velocity and other parameters which are not accurately 
known. The overall uncertainty in mass determination is estimated to be about 
an order of magnitude. The problem is then compounded by applying one of 
several empirical penetration formulas to the mass distribution to obtain a depth 
of penetration in a particular material. 

The present approach makes a more direct use of the observational data 
and does not require specification of mass for any observation, thus avoiding the 
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major source of uncertainty. It is necessary only to know the exponent of mass 
in the equations relating mass to the observed quantities. There is some 
theoretical justification for the values of these exponents, although in some cases 
empirical values have been found that differ slightly from the theoretical values. 
The method employed is a differential technique that makes use of two or more 
similar observations with different thresholds to obtain values for the logarithmic 
derivative d log N/d log m. It is important to distinguish between the observed 
number of incidents per unit area-time with mass, velocity, density, etc. , such 
that the observed quantity is above the threshold of observation, and the number 
N incident per unit area-time-solid angle with mass equal to or greater than 
some characteristic value. These quantities are related by a factor which is a 
function of d log N/d log m and contains integrations over angular, velocity, and 
density distributions. Fortunately, since present knowledge of such distributions 
is at best sketchy, they only enter in these correction terms, and therefore exert 
only a minor influence on the final result. 

Since it is desirable to avoid the necessity of determining mass for any 
observation, the logarithmic derivatives are specified at values of log N, which 
are taken as the mean of the logarithms of the fluxes that determine a particular 
log derivative. In this manner, a set of values for the log derivative and log N 
are found from the observations from which a functional relation F (log N) may 
be found by the method of least squares. The resulting differential equation 


d (log N) 
d (log m) 


F (log N) 


may be integrated to give 


( 1 ) 


log m 


r d log N 
J F (log N) 


+ const. 


( 2 ) 


By the previous assumption that the observable, in this case penetration 
of a thickness p of some specified material, is related to mass by 


log p = K + y log m (3) 

where K is an undetermined constant containing velocity, density, and material 
property terms and y is the known exponent of mass in the equation relating 
mass to the observable. Therefore, 


log p = 



d log N 
F (log N) 


+ const. 


( 4 ) 
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Assuming the velocity and density distributions are independent of mass, 
which appears to be at least approximately the case, the above equation holds 
over the range for which F (log N) is known. The constant of integration which 
contains all the uncertainties in the definition of mass, velocity, density, and 
material properties may be determined by a single experiment which yields a 
value of N for known p of a given material. The 0. 4-mm Pegasus datum point 
is used to obtain this constant and the derived frequency of penetration versus 
thickness of material applies to 2024-T3 A1 backed by a soft foam. 

The mass distribution function is found from the cumulative encounter 
frequency N. The differential equation derived from the measured log derivatives 
gives a relation of log m and log N to within a single constant. This constant is 
determined to give a best fit of the resulting values estimated from each meas- 
urement. This obviously does require the assumption of velocity and density 
distributions, and is also subject to errors in mass determination. However, 
since all the measurements are combined to find only one constant, hopefully 
some of these errors are averaged out, and the result is a mass distribution 
that is most consistent with the entire data set. 


FUNDAMENTAL RELATIONS 


Unfortunately, none of the present measurement techniques directly 
measure meteoroid mass. All such measurements have an observed quantity Q 
that is related to mass by an observational equation of the form 

log Q = log K + a log p + /? log v + y log m + 6 log cos 6 (5) 

where K is a constant, p the meteoroid density, v the meteoroid velocity, m the 
mass, and 9 the angle of incidence. For penetration experiments, Q is the 
thickness of the detector; for radar observations, Q is the electron trail density; 
and for photographic meteors, Q is the intensity or the photographic magnitude. 

Because of the spectrum of velocities, densities, and angles of incidence, 
there is no unique value of m corresponding to Q. The observed frequency for 
a given threshold Q is 


OO OO OO 
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where n (p) is the normalized density probability density function, 

P 

n (v) is the normalized velocity probability density function, 
v 

n (m) is the average isotropic directional mass flux differential distribution 
m 

function, and m (Q, p , v, 6) is the threshold mass having density p, 

velocity v, and incident angle 0 capable of producing Q. It is assumed that 

functions n (p) , n (v) , and n (m) are independent, 

p v m 

Thus, each observation may be written as an integral equation, the set 
of which may be solved for the mass distribution n^(m) . 


It is useful to define a characteristic mass m of a meteoroid having 

c 

normal incidence, average velocity <^v^> and average density ^p^> that is capable 
of producing a threshold Q. Obviously this is contained in the interval of m 

that are the major contributors to cp. Further, it may be assumed that this 
interval is small compared to the total range of m, primarily because the mass 
distribution decreases rapidly for larger masses, and the velocity and density 
distributions decrease rapidly at some upper bound preventing the larger pop- 
ulation of smaller masses from being observed. The intergral may be written 

in terms of N , the number of meteoroids with mass m or greater incident per 
m 

unit area-time -solid angle, 

oo 

N = f n ( m ) dm . (7) 

m J m 
m 


Expanding log in a Taylor T s Series about log m^, 


( d log N 

j-. — — 

d log m 
, /m 


[log m - log IMJJ+ 


Let 


'd log N 


X = 


m 


, d log m 


m 


( 8 ) 


( 9 ) 
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The minus sign is introduced for convenience so that A will be a positive 
quantity. 


Then 


f n (m) dm = 
J m 


m 


m 


m 


N 


m 


( 10 ) 


The observed flux (p may be written in terms of N the incidence flux per unit 

c 

solid angle of mass m or greater 

c 


N 


m 


</> = 


<p> “ Vr <V Vy 


ffcos 1 + 6x/r edafn (p)p“ x/r d P 


hemisphere 


/ V v > 

0 


v^dv . 


( 11 ) 


Integrating, 


<p = N 


m (l+(5A/2y) / 
c ^ 


<p aX/ ' l '> <v^> 

<P> aVr <vf X/y 


( 12 ) 


where m , from its definition, is 
c 


log m c = — (log Q - log K - a log^p^ - p log ^ ) . 


(13) 


The factors appearing in the relation between incident flux and observed 
flux may be combined into a correction term 


F (A) = 


7T 


(l + <5A/2y) 


(p aX/y ) 

( p ) ax/y < v >' JX/Y 


(14) 
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The first factor is an obliqueness correction which accounts for angles of incidence 
other than normal. The density and velocity terms adjust the observations ac- 
cording to the contributions from masses above and below the characteristic 

mass m . 
c 


Finally, the number incident having mass m 

c 

from the observed flux by 


or greater may be obtained 


N = <b/ F(A) . 
m 


The quantity A must be found from the log derivative 


A 


d log N 
m 

d log m 


d log <p d log F (A) 

d log m d log m 


(15) 


(16) 


This may be solved by the method of successive approximations using as 
a first approximation 


A o = d log 0 ^ _ A log 0 = log (flj/ift 2 ) 

d log m A log m 7 log (Qj/Q 2 ) ( ’ 

where <p± and <p 2 are measured rates which correspond to observables Qjand Q 2 
for which the same observational equation applies, and y is the mass exponent 
in the observational equation. 


DENSITY AND VELOCITY DISTRIBUTION 


The determination of a velocity distribution for meteors has been prob- 
lematic because of the strong dependence of the observable on velocity. This 
results in a distortion of the observed velocity distribution because of the con- 
tribution of the more numerous smaller meteors on the high velocity end of the 
spectrum. One method of eliminating this selection effect was proposed by 
Orrok in which observations are categorized into groups of constant mass de- 
termined from 

T 

m = f I dt (18) 

t„v 3 J 
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where To is the luminous efficiency, T is the lifetime, and I the measured in- 
tensity. Taking a group having sufficient mass to insure practically ail of the 
slower meteoroids will be observed should give a reasonably unbiased set of 
data from which a velocity distribution can be obtained. 


Dohnanyi [ 1] applied this technique to a group of meteor observations by 
McCroskey and Posen and obtains a distribution function expressed by 


n 

v 


(v) = 


(c v 1-S 

| C x i. 61 x 10 7 v" 4 - 3 , 


11. 2 £ v ^ 16. 6 km/sec 
16. 6 < v < 72. 2 km/ sec 


(19) 


where C is the normalization constant. This distribution results in <^> = 19. 17 
km/sec. Figure 1 shows this representation. 


For lack of more definitive data or contradictory evidence this dis- 
tribution will be assumed to apply to the entire mass spectrum of meteors. 

The density distribution is more difficult to determine. Generally, the 
most definitive data available is the average density of a group of observations. 
Jacchia, Verniani and Briggs [2] find an average p of 0. 26 gm/cm 3 for photo- 
graphic meteors. Verniani and Hawkins [ 3] for smaller masses find relative 
abundances between several groups of densities. Their tentative findings are 
51 percent in a group less than 1 gm/cm 3 , with an average of 0. 37 gm/cm 3 , 45 
percent between 1 gm/cm 3 , with an average of 2. 8 gm/cm 3 , and 4 percent 
greater than 12 gm/cm 3 . The latter densities are of course unrealistically high, 
and probably result from small inaccuracies in the observed decelerations. For 
the lack of more definitive data, the assumed density distribution will be taken 
as 


n (p) = 0. 51 x 0. 37 6 (p - 0. 37) + 0. 45 x 2. 8 5 (p - 2. 8) 


+ 0. 04 x 8 <5(p - 8) 


( 20 ) 


where 6 (p - pj) are Dirac-Delta functions. The value of 8 was arbitrarily 
assigned to the high density component reported by Verniani which represents 
the highest density which could be reasonably expected. Since the observables 
dealt with are fairly insensitive to density and density distribution, the inac- 
curacies in this distribution do not significantly affect the results. This choice 
of distribution gives (py = 1. 77 gm/cm 3 which is somewhat higher than usually 
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assumed. This is due to the higher density contributions which may indeed be 
significant for the smaller particles. Figure 2 shows the assumed distribution. 


OBSERVATIONAL DATA 
Satellite Penetration Measurements 


For determining the characteristic mass for penetration of a given de- 
tector, the empirical Manned Space Center penetration formula has been utilized. 
This formula is usually written for semi-infinite targets in the nondimensional 
form 


= 1. 64 d 1//l8 (cm) (Pp/Prp) 1 ^ 2 (v cosa/c) 2 / 3 (21) 

where p is depth of penetration 

d is equivalent projectile diameter 
Pp is projectile density 

P T is target density 

v is projectile velocity 
c is target sound velocity 
9 is the angle of incidence 

The constant 1. 64 applies to relatively hard (2024-T3) aluminum. It is 
more convenient to express this formula for aluminum in the form 

p(cm) = 0.41 m°- 352 (gm) p°' 147 (gm/cm 3 ) v 0-667 ( 22) 

(km/sec) cos 2 / 3 9. 
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v (km/sec) 

FIGURE 1. THE VELOCITY PROBABILITY DENSITY 
FUNCTION OBTAINED BY DOHNANYI 



^(gm/cm 5 ) 


FIGURE 2. THE RELATIVE ABUNDANCES ASSUMED 
FOR THE VARIOUS REPRESENTATIVE DENSITIES 



To correct for finite target thickness, the constant 0. 41 is increased by a factor 
of 1.5. Thus the penetration formula may be put in the specified form for the 
observational equation by assigning 

Q = depth of penetration (cm) 

K = 0. 605 for 2024-T3 A1 

a = 0 . 147 

p = 0.667 

y = 0.352 

6 = 0.667 

For the Cu-Be and stainless steel pressure cells used on Explorer XVI 
and Explorer XXIII, a factor of 2 was used to correct for finite thickness be- 
cause the material is under hoop stress, has no backing, and requires only a 
leak for detection. Laboratory tests at AEDC [ 4] in stainless steel indicate a 
factor of 1. 47 less penetration in stainless steel than in A1 for velocities up to 
8 km/sec. Calculations by Bjork [5] indicated a similar value for velocities of 
20 and 72 km/sec. Tests at Langley Research Center indicate a similar re- 
sponse in Cu-Be and stainless steel at velocities up to 6 km/sec. Therefore, K 
is taken as 0. 41 x 2/1. 47 = 0. 558 for both Cu-Be and Stainless Steel pressure 
cells. 


E xplorer XVT . - During the useful life of Explorer XVI, 44 punctures were 
observed in the 0. 0025-cm Cu-Be pressure cells and 11 events in the 0. 0050-cm 
cells. The rates are 0. 333/m 2 day and 0. 172/m 2 day, respectively [6] . The 
semimajor axis is 965 km, which results in an earth shielding factor of 0. 748. 
The unshielded puncture rates are 0. 445/m 2 day and 0. 230/m 2 day, respectively. 

From these data: 


A o = _ y log (<ft i/<f>2) 
log (Q1/Q2) 

J = n /^7 = 10- 0 - 495 

F(X°) = 2. 21 
N = 0/F(A°) = 10"°- 839 


0. 334 ± 0. 15 
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Explorer XXIII. - In one year Explorer XXU3 observed 49 punctures of 
the 0. 0025-cm stainless steel pressure cells and 74 punctures in the 0. 0050-cm 
cells. The rates are 0. 345/m 2 day and 0. 205/m 2 day, respectively [7] . The 
semimajor axis of the orbit is 390 km, which results in an earth shielding 
factor of 0. 674. Removing this results in puncture rates of 0. 526/m 2 day and 
0. 305/m 2 day. 

From these data: 

A. 0 = 0. 278 ± 0. 11 

= io -0 - 397 

F(A°) = 2.32 

N = io ” 0 - 763 

Pegasus. - Since the three Pegasus satellites had similar orbits, their 
data were combined. As of January 1, 1966, the 0. 0038-cm detectors had 
counted 582 events for a rate of 0. 188/m 2 day; the 0. 02-cm detectors had counted 
49 events for a rate of 0. 021/m 2 day; and the 0. 04-cm detectors had counted 201 
events for a rate of 0. 00487/m 2 day. 1 The shielding factor is 0.69. The un- 
shielded rates are 0. 273, 0. 0305, and 0. 00706/m 2 day, respectively. The 
aluminum thicknesses quoted are front sheets of parallel plate capacitors which 
are temporarily shorted by meteoroid penetration. In addition, a 0. 0012-cm 
mylar dielectric must also be penetrated. The effect of this mylar in terms of 
equivalent thickness of A1 is not yet known, but since it is thin compared with 
the 0. 02 and 0. 04 plate, no correction is applied to these thicknesses. The 
mylar is a significant fraction of the 0. 0038-cm detector and must be considered 
in this case. Also, the 0. 0038-cm detector is 1100-6 aluminum rather than 
2024-T3, and is backed with a 0. 013-cm hard epoxy layer rather than soft foam. 
An effective thickness in terms of the 0. 02- and 0. 04-cm detectors has not yet 
been accurately determined; therefore, this datum point will not be used to 
determine the slope of the penetration curve. From the 0. 02-cm and 0. 04-cm 
points , 


A 0 = 0. 740 ± 0. 14 

J = io - 1 - 841 


1. Clifton, K. S. and Naumann, R. J. : Pegasus Satellite Measurement of 
Meteoroid Penetration. NASA TM X (in publication). 
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1. 76 


F(X°) = 

N = 10~ 2 - 076 


Radar Observations 


Elford [ 8] using a transmitted power of 1. 2 megawatts, reported a total 
flux of 9. 05 x 10 -4 /m 2 day of meteors whose maximum line densities exceed 
3. 3 x 10 10 e/ m. Reducing the transmitted power by a factor of 143 resulted in a 
decrease in observed echoes by a factor of 12. Assuming the reflected electric 
field vector is proportional to the average trail density, the transmitted power 
required for a given reflected power varies as the inverse square of the electron 
line density. Thus, the minimum detectable line density for the low power 
operation is 3. 94 x 10 11 e/m. 

Using the ionization efficiency determined by Verniani [ 9] , the maximum 
line density (e/m) may be related to mass by 

log q = 10.638 + log m + 4 log v + cos 0 (23) 

where v is in km/ sec and m is in grams. 

The observational equation may be specified by setting 

Q = electron line density, q (e/m) 

log k = 10. 638 

a =0 

13 =4 

y =1 

6 =1 

The earth shielding factor is unity since the area in question is parallel 
to and facing away from the earth. 
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The first approximations are 


A. 0 = 1. 005 ± ? 

0 = i 0-3. 583 

F(A°) = 5. 94 

N = 10“ 4 - 358 


Photographic Observations 


Hawkins and Upton [10] determined the influx rate of photographic meteors 
in the range of -2. 5 to +3 photographic magnitude. Their observed flux may be 
expressed in terms of peak photographic magnitude Mp or brighter as 

log <p = 0. 537 Mp - 8. 95 (no./m 2 day) . (24) 

The uncertainty in this relation is least from -2. 0 < Mp s 2. 5. 

Jacchia, Verniani, and Briggs [ 3] performed a detailed analysis of a 
selected group of Super-Schmidt meteors whose masses were determined by 
integration of the total light curve. They found by method of least-squares a 
relation for peak brightness given by 

Mp = 11.59 - 8. 75 log v - 2. 25 log m - 1. 5 log Q . (25) 

log Q = Mp 

log K = 11.59 

a - 0 

/3 = -8.75 

y = -2. 25 

8 = - 1.5 
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From these data 


A° = 1. 205 ± 0. 06 

^ = 10 - 8-815 

F(A°) = 10.41 

N = io -9 - 832 


COMPUTATION PROCEDURE 


The correction factor F(A) which relates observed number to cumulative 
directional flux depends on the measurement involved. This factor is shown in 
Figure 3 for penetration, radar, and photographic measurements. A few comments 
concerning the interpretation of this factor may be in order. The factor con- 
tains 7 r which, when multiplied by the number per steradian, gives the total 
number incident on a flat surface. For penetration observations, the fact that 
F(A) is less than 7 r means that most of the observed punctures are caused by 
meteoroids with mass greater than the characteristic mass. This is primarily 
because of the reduced penetrating ability of particles impacting obliquely. For 
the radar and photographic observations, the large velocity dependence allows 
a significant number of the more numerous meteors with masses smaller than 
but with velocities higher than average to be observed; hence, the value of 

F(A) is greater than 7 r . 

The first approximations for A versus y = -log N are shown in Figure 4. 
Successive approximations may be obtained by the following process. Taking the 
log of equation (15) and differentiating, 

d log N = d log cf> - ^° gF(X) ~ n d log N (26) 

Let x = log m, y = -log N 

d£ = d log eft + d log F (A) dA dy_ 

dx dx dA dy dx 
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FIGURE 3. THE RATIO F(X) OF THE OBSERVED FLUX 

OF METEOROIDS TO THE NUMBER ENCOUNTERED PER FIGURE 4. THE POSSIBLE FITS TO X AS A 

UNIT AREA TIME-SOLID-ANGLE WITH MASS SUF- FUNCTION OF Y. THE DEPARTURES OF THE 

FICIENT TO PRODUCE A THRESHOLD OBSERVATION ITERATED VALUES OF FIT I FROM THE OB- 

PROVIDED THEY HAVE AVERAGE DENSITY, AVERAGE SERVED VALUES ARISE FROM THE COR- 

VELOCITY, AND NORMAL INCIDENCE RECTIONS IN X INVOLVING dX/dy 




By definition, 


A = -*• ; A» = , x , a»/ 

dx dx 


[" _ d log F(A) dxl 
clA dyj • 


(28) 


A functional relation for A 1 = A(y*) is found by the method of least squares 


and is used to compute dAVdy 1 . The (i + l) u ** approximation is found for each 
measurement from 


th 


i + 1 


th 


= A 0 / 


1 - 


d log F1A 1 ) dA 1 


dA 


dy 1 -I 


(29) 


The (i + 1) approximation for each y is then found from 


i + 1 i + 1 

y = log F (A ) - log 0 


(30) 


i + 1 i 

The process continues until A = A for each value of A. 


The question of the analytical representation for A = A(y) must be con- 
sidered. The most straightforward choice that could provide a reasonable fit 
is a second degree polynominal, 

A = ay 2 + by + c (31) 

It may be seen from Figure 3 that d log F(A)/dA is negative for penetration ob- 
servations for A < 1.1, and that it is positive for radar and photographic obser- 
vations in the range of A applicable. Since the first approximations for A in 
Figure 4 indicate positive but decreasing value of dA/dy for increasing y, it may 
be seen from equation (28) that the values of A for penetration observations will 
be decreased while those corresponding to radar and photographic observations 
increase. Because of the larger dA/dy at the radar point, the corresponding A 
is increased more than the A corresponding to the photographic point. This 
causes the polynomial fit to reach a maximum before the photographic point and 
dA/dy becomes negative at the photographic point. This change of sign forces 
the corrected A for the photographic point to have a smaller value, and the final 
result is that indicated by Fit I in Figure 4. Physically, this may be interpreted 
to mean that the mass distribution is such that the higher velocity tail of the 
velocity distribution, which allows the more numerous smaller mass meteors to 
be above the observational threshold, is less effective for the larger masses 
than for the smaller masses. Hence the observed slope d log <p/d log m is 
steeper than the true slope d log N/d log m. 
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There is, of course, no theoretical justification in the choice of the 
polynomial chosen to represent X(y) . The peak in X between the radar and 
photographic data may not correspond to physical reality. For this reason, an 
alternative representation was chosen that conditions X(y) to be monotonic, i. e. , 


i 

X 


a + 


y - c 


(32) 


The results of this choice are indicated as Fit I in Figure 4. There appears to 
be considerable difference in the two results, which more or less represent 
extremes in the possible choices for A(y). Fortunately, however, the end 
results, such as finding a mass distribution or penetration frequency, are not 
seriously influenced by the choice of representation of X(y) . 


The cumulative mass distribution may be found to within a constant of 
integration by solving the differential equation 


dx 


= My) 


(33) 


or 

logm = / ^ + d (34) 

For Fit I, 

_2 

log m = . — tanh -1 

'v/b 2 - 4ac 

where from the fit of X = A*(y) 

a = -0. 0328156 
b = 0. 425222 

c = -0. 550693 

For Fit H, 

log M = - a log N - b In (-log N-c) + d^. (36) 


-2 a log N + b 
N/b 2 - 4 ac 


(35) 
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where a, b, c have been determined from the fit of X = A^(y) 

a = 0.719683 
b = 0.858195 
c = 0. 590827 

A value for d was computed for each datum point. Table I is a summary 
of these results. The average value and RMS deviation in d is -2. 35 ± 0. 53 for 
Fit I and -8. 45 ± 0. 46 for Fit II, with the radar points producing the largest 
deviations in both cases. Using the average value for d, the cumulative mass 
distribution is shown in Figure 5. It may be seen that the satellite points fall 
slightly to the right of the mass distribution and photographic points fall slightly 
to the left. This means that meteoroids are either somewhat more penetrating 
than predicted from the MSC formula or that the masses inferred from photo- 
graphic data are too small. The radar points fall considerably on the small 
side. However, the uncertainty of ± 0. 3 in the velocity exponent of the radar 
observational equation estimated by Verniani [ 3] is sufficient to bring these 
points into agreement. 

It is not necessary to actually use the mass estimates with their inherent 
uncertainties to obtain an expression for penetration frequency. Since penetration 


P ~ m 7 , (37) 

log p = y log in + constant (38) 

log p = -1. 6897 tanh -1 [ 0. 1575 log N + 1.0206] + d* (39) 

for Fit I, or for Fit II, 

log p = -0. 2533 log N - 0. 3021 In (-log N - 0. 5910) + d* (40) 


The constant d* now contains the combination of material velocity, and 
density terms in the penetration formula as well as the constant of integration in 
the mass distribution. This constant can be determined for a particular material 
from a single penetration measurement. Using the Pegasus 0. 04-cm datum 
point, 

d* = -0. 1054 
d* = -2.060 
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TABLE la 


Fit I 

Measurement 

<p Obs. 
no./m 2 day 

A 

F(A) 

N 

no./m 2 day sterad 

log d 

Est. Mass 
(gi») 

Mass From Fit 
(g™) 

.04 cm. Peg. 

i 0 -2. 16 

.773 

10-230 

Iq-2.39 

-2.230 

10-5.91 

Iq-6.03 

. 02 cm. Peg. 

1 0-1. 53 

. 601 

10-264 

10-1.79 

-2.211 

10-6.76 

10-6. 90 , 

. 0050 cm. Exp. 
XXHI 

10" 516 

. 292 

Iq.360 

1Q--876 

-1.673 

10-8.33 

10" 9 - 01 

. 0025 cm. Exp. 

xxni 

iO-27S 

. 215 

1q. 392 

10". 670 

-1.709 

10-9.19 

10-9.83 

. 0050 cm. Exp. 
XVI 

10 -.638 

.331 

10-345 

10“ -983 

-2.017 

Iq-8.33 

10-8-67 

. 0025 cm. Exp. 
XVI 

1Q-.352 

.239 

10-382 

10". 734 

-1.990 

10-9.19 

Iq-9.55 

Low Power Radar 

i0-4.12 

1. 270 

IqI.092 

10-5.21 

-3.209 

10-4-18 

CM 

CO 

CO 

1 

o 

t-H 

High Power Radar 

10 -3,04 

1. 109 

10-888 

^ q “3 # ^ 

-3.210 

10~5.25 

IQ -4 * 40 

-2 Photo. Mag. 

IQ-10. 02 

. 760 

10-597 

i 0-10.62 

-2.686 

10 1 * 05 

10 1 ' 38 

2. 5 Photo. Mag. 

10-7-61 

1. 176 

10-98! 

Iq-8.59 

-2.608 

10--948 

1Q-.694 




Fit II 


Measurement 


< p Obs. 
no./m 2 day 



TABLE lb 


no./m 2 day sterad log d 


, Est. Mass 


Mass From Fit 
(gm) 








FIGURE 5. THE CUMULATIVE MASS FLUX DISTRIBUTION THAT BEST 
FITS THE ESTIMATED CHARACTERISTIC MASS FOR EACH MEASUREMENT. 
THE ERROR BARS ON THE RADAR POINTS IN THE VICINITY OF 10“ 5 
gm INDICATE THE RESULT OF UNCERTAINTY OF THE DEPENDENCE 
OF IONIZATION EFFICIENCY ON VELOCITY 
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PENETRATION FREQUENCY 


Figure 6 is a plot of the estimated frequency of penetration versus thick- 
ness of 2024-T3 A1 backed with foam. These same curves may be used for any 
other material by establishing an equivalent thickness of 2024-T3 for the material 
by laboratory impact tests. The pessimistic curve results from Fit II. The 
maximum difference between Fit I and II is a factor of 0. 4-in. thickness re- 
quired to give a certain penetration frequency. It is recommended that the largest 
value be used in all cases for a conservative design value. Table H is a summary 
of penetration frequency for various thickness of 2024-T3 Al. 

The other possible source of uncertainty in the penetration formula lies 
in the values of y used in obtaining the log derivatives. Recall that empirical 
values of y = 0. 352 for penetration calculations and y = -2. 25 were used for 
photographic analysis. These are slightly different from their respective 
theoretical counterparts of 0. 333 and -2. 5. 

To investigate the possible effects of this uncertainty, the calculations 
were also run with these values. The optimistic curve results from Fit I with 
the smaller values of y, i. e. , 0. 333 and -2. 5. Since the two major sources of 
uncertainty, i. e. , the choice of analytical representation and the mass exponents 
have been analyzed with more or less extreme cases, the envelope of the re- 
sulting penetration curves may be considered to represent the limits with which 
penetrating flux can be predicted with available data. Although the limits diverge 
with increasing thickness , the spread is only a factor of 2 in thickness for the 
largest thickness of interest. Other choices of velocity or density distributions 
have virtually no effect on this result. 

For thicknesses less than 0. 04 cm, the difference between the various 
cases is less. Figure 7 is an expanded plot of this region. The problem of 
material equivalence between the 0. 0038-cm Pegasus point, the Explorer XVI 
points, and the Explorer XXIII points is avoided by the fact that the technique 
uses the slopes and the number observed rather than the actual thickness. Values 
for material equivalence can be inferred from this plot by taking the thickness in 
2024-T3 Al from the plot that corresponds to the observed penetration frequency. 
The results are: 
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FIGURE 6. ESTIMATED PENETRATION FREQUENCY IN VARIOUS 
THICKNESSES OF 2024-T3 AI. THE ENVELOPE REPRESENTS LIMITS 
OF COMBINED UNCERTAINTIES IN MASS DEPENDENCE AND FORM OF 
FUNCTIONAL FITS. THE INDICATED MEASUREMENTS ARE PENETRATION 
RATES ESTIMATED ON THE BASIS OF THE CALCULATED CHARACTERISTIC 

MASS FOR EACH MEASUREMENT 


t . 
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TABLE II 


Thickness 

Upper Limit 

Lower Limit 

>4_T3 Plate (cm) 

no. /m 2 day 

no. /m 2 day 

. 003 

4.40 x 10 -1 

3.60 x 10 _1 

. 004 

3. 52 x 10 _1 

2. 72 x 10 _1 

. 005 

2. 81 x 10 -1 

2. 15 x 10' 1 

. 006 

2. 31 x 10 -1 

1. 70 x 10 -1 

. 007 

1. 92 x 10" 1 

1. 40 x 10 -1 

. 008 

1. 63 x 10- 1 

1. 19 x 10 _1 

. 009 

1. 49 x 10 -1 

1. 01 x 10 -1 

. 010 

1. 19 x 10 -1 

8. 90 x 10 -2 

. 020 

3. 35 x 10 -2 

2. 90 x 10 -2 

. 030 

1. 39 x 10" 2 

1.28 x 10 -2 

. 040 

7. 10 x 10" 3 

7. 10 x 10~ 3 

. 050 

4. 30 x 10“ 3 

4. 30 x 10 -3 

. 060 

2. 50 x 10 -3 

2. 36 x 10 -3 

. 070 

1. 77 x 10~ 3 

1. 60 x 10“ 3 

. 080 

1. 28 x 10“ 3 

1. 10 x 10 -3 

. 090 

8. 50 x 10 -4 

7. 40 x 10' 4 

. 100 

6. 50 x 10~ 4 

5. 40 x 10“ 4 

. 200 

8. 50 x 10“ 5 

5. 60 x 10~ 5 
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TABLE II (Concluded) 


Thickness 
2024-T3 Plate (cm) 

Upper Limit 
no./m 2 day 

Lower Limit 
no. /m 2 day 

. 300 

2. 50 x 10 -5 

1. 20 x 10 -5 

.400 

9. 50 x 10“ 6 

4. 00 x IO -6 

. 500 

4. 98 x 10" 6 

1. 65 x 10 -6 

. 600 

2. 70 x 10 -6 

7. 90 x 10 -7 

. 700 

1. 61 x 10 -6 

4. 10 x 10“ 7 

. 800 

1. 09 x IO -6 

2.40 x IO" 7 

. 900 

7. 30 x 10 -7 

1.46 x 10“ 7 

1. 000 

5. 20 x 10 -7 

9. 50 x 10" 8 

2. 000 

5. 40 x IO" 8 

5. 80 x 10" 9 

3. 000 

1. 30 x IO" 8 

1. 15 x 10~ 9 

4. 000 

5. 02 x 10“ 9 

4. 00 x IO -10 

5. 000 

2. 40 x IO" 9 

1. 89 x IO -10 

6. 000 

1. 29 x 10 -9 

1. 10 x IO" 10 

7. 000 

7. 40 x 10 -1 ° 

6. 80 x IO -11 

8. 000 

4. 60 x 10 _1 ° 

4. 50 x IO -11 

9. 000 

3. 00 x 10" 10 

3. 10 x IO -11 

10. 000 

2. 10 x IO -10 

2. 30 x IO -11 



Detector Material 
and thickness (cm) 

Equivalent thickness 
in 2024-T3 Al (cm) 

Equivalence 

Factor 

0. 0025 Cu-Be 

0. 0028 ± 0. 0009 

1. 12 ± 0. 36 

0. 0050 Cu-Be 

0. 0056 ± 0. 0020 

1. 12 ± 0. 40 

0. 0025 Stainless 

0. 0021 ± 0. 0009 

0. 84 ± 0. 36 

0. 0050 Stainless 

0. 0042 ± 0. 0010 

0. 84 ± 0. 20 

0. 0038 1100 Al 

0. 0046 ± 0. 0007 

1. 21 ± 0. 18 


It may be seen that all the sensors have material equivalence factors near unity, 
which means they are essentially equivalent to 2024-T3 Al. It is known that 
both stainless steel and Cu-Be semi-infinite target are penetrated less deeply 
by hypervelocity projectiles. However, apparently this is compensated some- 
what by the fact that these detectors have a free rear surface which is more 
subject to spall than a foam-backed material. Recall that the estimated equiv- 
alence factor (Page 11) was 1. 11, which is in good agreement with the derived 
value for Cu-Be. It appears that the stainless steel is not as puncture resistant 
as the Cu-Be since the puncture rates in both thicknesses of stainless steel 
exceed the corresponding rates in Cu-Be, and the slope determined from the 
data points is not as steep for Explorer XXIII. Of course the small statistical 
sample does not permit this to be stated definitely. The difference between the 
observed rates in the 0. 0025-cm detectors is 18 percent while the standard 
deviations associated with the rates are approximately 14 percent. Also, the 
standard deviations in the slope are not sufficiently small to say that the slope 
corresponding to stainless steel is actually less than that of the Cu-Be. The 
0. 0038-cm Pegasus detector has a 0. 0012-cm mylar backing that also must be 
penetrated to produce a signal. Thus, the total material thickness is actually 
0. 005 cm which is quite close to the 0. 0046-cm equivalent thickness. Apparently 
the fact that 1100 Al has lower strength than the 2024-T3 is compensated by the 
thick epoxy backing on the 0. 0038 cm-Pegasus detector. 


MASS DISTRIBUTION 


As was seen in Figure 5, the cumulative mass distribution obtained from 
either fit agrees with the values estimated from the various measurements to 
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within ± approximately half an order of magnitude. Both representations produce 
computed values of A that agree with the observed values to within observational 
limitations. 



A 0 OBS 

A 0 Fit I 

A 0 Fit n 

Explorer XVI 

0. 334 ± 0. 15 

0. 324 

0. 331 

Explorer XXIII 

0. 278 ± 0. 11 

0. 290 

0. 264 

Pegasus 

0. 740 ± 0. 14 

0. 724 

0. 794 

Radar 

1. 006 ± ? 

1. 007 

0. 980 

Photographic 

1. 205 ± 0. 06 

1. 160 

1. 206 


The Fit I gives slightly less difference between the constant of integration 
determined from the Pegasus points and that determined from the photographic 
points. This difference is 0. 47 in the logarithim or a factor of 3 in mass. It is 
reasonable to assume that such a difference can be accounted for by our imperfect 
knowledge of penetration mechanics, particularly for threshold penetrations in 
finite targets, and by the uncertainty inherent in inferring mass from photo- 
graphic observations. 


The differential directional mass flux distribution function n (m) can be 

m 

constructed for any value of mass m^ from 


n (m ) dm = N , - N 

m o m - dm m 

o o 


(41) 


using the Taylor expansion for about m^, equation (10) , 


n (m ) dm = 
m ' o 


(m - dm 

o 

V m 


-A 


N - N 
m 


m 


(42) 


n (m ) = 
m ' o 


N 

m m 


(43) 
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The resultant mass distribution is shown in Figure 8. It is interesting 
to compare this derived mass distribution with that of Van de Hulst [ 11] obtained 
from interpretation of the F-component of the solar corona and zodiacal light in 
terms of light diffracted from dust particles in interplanetary space. Van de 
Hulst finds the photometric data can be well represented by a distribution in 
radii by 


N(a) = c a -2,6 (no./ cm 3 ) 


(44) 


where 


c = 


10 - 2 ° 

3. 5 x 10" 19 
0 


at 1 A. U. 
at 0. 5 A. U. 
<0. 1 A. U. 


for a < 0. 035 cm. 


Taking the average density as 1. 7 gm/cm 3 and the average velocity as 
20 km/ sec, this distribution can be transformed to a distribution in directional 
mass flux encountered by an object in an orbit similar to that of the earth, 


n^fm) = 1. 25 x 10 6 m -1 - 534 (No./m 2 day str. gm). 


• It may be seen from Figure 6 that the agreement both in magnitude and 
slope between the mass distributions is remarkable for masses from 10 -5 gm to 
10 _1 ° gm. Van de Hulst states that particles with radii < 0. 035 cm are less 
abundant than indicated by his distribution function, which agrees very well with 
present result. 

A flattening and possible eventual decline of the distribution function for 
smaller masses would be expected as the result of the increased significance of 
radiation perturbations such as Poynting-Robertson effect and radiation pres- 
sure. The derived expressions for N approach a maximum value for the incident 
number which is 0. 738/ m 2 day steradian for Fit I and 0. 255/m 2 day steradian for 
Fit H. The existing data do not extend to small enough masses to determine 
this accurately, but it is reasonably certain that only approximately 2-3 en- 
counters /m 2 day can be expected with all meteoroids regardless of their mass. 
This is consistent with the fact that only one encounter in 1. 4/m 2 day exposure 
has been observed on Gemini windows 4 which sets a 90 percent upper limit of 
2. 8/m 2 day on events that produce discernable damage to optical surfaces. It is 

i. Burbank, P. : Manned Space Center, private communication, 1966. 
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FIGURE 7. AN EXPANDED REGION OF THE PENETRATION 
CURVE. THE EXPLORER POINTS AND THE 0. 038 mm PEGASUS 
POINT WERE ADJUSTED IN THICKNESS TO AGREE WITH THE 
SLOPE OF THE PENETRATION CURVE. THE RESULTING 
VALUE IS TAKEN AS THE EQUIVALENT THICKNESS IN 

2024-T3 AI 


FIGURE 8. THE DIRECTIONAL MASS 
FLUX PROBABILITY DENSITY FUNCTION 
COMPARED WITH ESTIMATES OBTAINED 
FROM ZODIACAL LIGHT BY Van de HULST 




also interesting to note that meteoroid measurements from Mariner IV [ 12] 
from 1 - 1. 56 A. U. observe between 0. 63 and 2. 8 events /m 2 day. 

No value for cut-off mass due to radiation pressure if one indeed exists 1 
can be inferred from the functions derived from either Fit I or II because the 
mass corresponding to the maximum N is 0. The mass contribution m n^ 

(m) dm from the particles smaller than 10“* 3 gm is negligible. The total 
mass influx may be found by integrating 

oo 

m_ = ff cos 6 dfi | m n (m) dm . (45) 

i X J m 

2 ir 0 

XN 

, m 

Since n (m) = 

m m 

OO 

rn = 7r f X N dm . 

T i 


For a small interval (mj, m 2 ) the contribution in mass influx is found approx- 
imately using the Taylor expansion for about m t 


m 12 



“l 



m 12 


7r X N m^ 
m l 1 
1 - X 



(46) 


(47) 


The total mass influx is found by summing the incremental contributions. 
This presents no problem for Fit II since X approaches an asymptotic value 
greater than unity. However, using Fit I, X approaches 0 for masses greater 
than those corresponding to the photographic data. This leads to an infinite 
mass influx if the integration is carried to m = ». Since such masses exceed 
the limits of the data, there is no justification for extending the derived expression 


1. Often this limit is estimated by equating gravitational force to radiation 

pressure applied to the geometrical cross section of a sphere and solving for 
the radius. This does not hold for dielectric particles or particles whose 

dimensions are of the order of or less than a visible wavelength. In fact 
Shapiro showed that submicron metallic andallsizesof Si0 2 particles never 
have radiation pressure forces greater than gravitational forces. 
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for N into this region. Therefore the summation is carried only to the last 
datum point. The mass influx contributions for larger masses are believed to 
be insignificant because of the rapidly decreasing population. If, for example, 
the value for A. for masses greater than 10 gm is 1. 2, the contribution from 
larger masses is only 3 percent of the total. 

The total mass influx is found by 2. 7 x 10 -7 gm/m 2 day for Fit I and 
1. 4 x 10 -7 gm/m 2 day for Fit II. This leads to a mass accretion of around 100 
metric tons per day. This is considerably less than earlier estimates based on 
extrapolation of photographic meteors to smaller masses which ranged from 
1000 to 10 000 tons per day. 


Having determined the mass distribution, there are several other quan- 
tities of interest that may be easily found. The normal component of momentum 
flux may be found from 


mv = //dfi cos 6 I dv v cos 0n (v) f dm m n (m) 
1 2n 0 V 0 m 


(48) 


mv = — rn 
_L 3 T 


<*> 


mv^= 0. 3 dyne - sec/m 2 day 

or a pressure of 2 x 10 -11 N/m 2 which is 10 5 times smaller than radiation pres- 
sure. 


The energy flux is 



E = 4 x 10 2 Joules/ m 2 day 

which over the entire earth amounts to a power input of 200 megawatts. 


COMPARISONS WITH OTHER EXPERIMENTS 


It appears that satellite penetration measurements together with radar and 
photographic observations of meteors yield a consistent picture of the earth's 
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meteoroid environment. There is considerable disagreement between this 
picture and the previously accepted meteoroid environment based on collection 
experiments and acoustical and photomultiplier sensing of impacts from very 
small meteoroids. This disagreement had been suspected for some time be- 
cause of the failure of various early experiments using wire grids or similar 
sensors that required physical damage to detect anywhere near the number of 
events predicted by the acoustical measurements. The low puncture rate ob- 
served by Explorer XVI confirmed this disparity. 

As an example of the order of magnitude of this discrepancy, the mass 
flux relation obtained by Alexander et al [ 13] from acoustic measurement of the 
momentum imparted to sounding boards is 

log N = -17 - 1. 7 log M 

where N is the number of encounters per square meter second with meteoroids 
of mass M or greater. The characteristic mass for the Explorer XVI 0. 0025- 
cm pressure cells is 10 -9 * 55 gm. From Alexanders results the cells should be 
punctured at a rate of 10 -1 * 5 per m 2 second or 10 3 * 4 per square meter day. This 
is almost four orders of magnitude higher than what was actually observed. 
Furthermore, the quantity d log N/d log M is -1. 7 according to Alexanders re- 
sults but was observed to be only -0. 3 from the penetration measurements. It 
is inconceivable that a measurement as straightforward as the Explorer XVI 
pressure cell experiment could measure four orders of magnitude too low. Most 
of the cells eventually registered a puncture, which is convincing evidence that 
the cells retained pressure and therefore were not punctured prior to the indicated 
time. There appears to be only two possible explanations of this discrepancy; 

(1) the acoustical data are in serious error, or (2) meteoroids are far less 
penetrating than expected. Although the empirical penetration formula used to 
predict characteristic mass for puncture is an extrapolation, it is based on 
laboratory data which contains the size range in question and velocities well in- 
to the meteoroid range (> 10 km/sec) . Most of the small mass, high velocity 
penetration data have been restricted to densities of 1 gm/cm 3 and higher, but 
theoretical studies [5] indicate high velocity particles of porous Al (0. 44 gm/cm 3 ) 
are just as penetrating for a given mass as those of normal density. If the 
acoustical measurements are actually measuring momentum of impacting mete- 
oroids, the only explanation of the failure to observe a large number of penetra- 
tions in very thin metallic foils must be a combination of very low relative ve- 
locities and extremely tenuous particles. Even if such particles were in geocentric 
orbits similar to the orbit of the satellite carrying the measurement, the fact that 
the satellites involved have fairly high inclinations requires the particles to have 
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relative velocities as high as 8 km/sec, with an average velocity of around 5 
km/ sec. The low penetration rate observed on the 0. 0025-cm Explorer XVI and 
xxni pressure cells and the failure to observe any penetrations on Jennison's 
[4] 0. 0012- and 0. 00l5-cm foils on Ariel II require incredibly low densities. 

It is clear that to explain both the acoustic measurements and the penetration 
measurements, it is necessary to assume two distinct families of meteoroids. 

The mass distribution obtained in this work applies to moderate density, 
penetrating particles in heliocentric orbits. The slope and magnitude of this 
distribution is such that a very small fraction of the acoustic events can be 
attributed to this family. If the acoustic events are actually caused by meteoroid 
encounters , the particles must belong to a family of extremely tenuous , non- 
penetrating particles in a distribution whose population declines rapidly with 
increasing mass. The much smaller decrease in population with mass observed 
for the family of penetrating particles indicates that an insignificant fraction 
of observed penetrations can be attributed to the more massive members of 
the family of non-penetrating particles. 

The assumption of this near-earth dust belt whose mass distribution is 
approximately four orders of magnitude higher than the interplanetary back- 
ground raises a major theoretical objection which has not been resolved despite 
considerable efforts on the part of several investigators. Because of their small 
size and apparent low density, such particles in a near-earth orbit could not 
remain indefinitely because of radiation perturbations. Therefore, they must 
somehow be supplied. An extensive analysis of possible concentration mechanisms 
was performed in which no means of providing a concentration of dust in the 
vicinity of the earth to the extent implied by the acoustic measurements was 
found. 1 


The difficulties encountered in attempting to explain the existence of such 
a near-earth concentration of dust as well as the peculiar properties that must 
be ascribed to such dust particles to account for their failure to penetrate very 
thin metallic foils demand that the evidence for their existence be reevaluated. 

The bulk of this evidence has been supplied by acoustical measurements. It has 
recently been demonstrated that certain types of acoustical impact measuring 
devices are susceptible to noise during thermal cycling. 2> 3 This could account for 
the fact that acoustical detectors operating in the vicinity of earth where they 
traverse the earth's shadow exhibit a high count rate while the same detectors 
operating in a fairly constant thermal environment on space probes give results 

1. Shapiro, 1. 1. ; Lautman, D.A. ; and Colombo, G. : The Earth's Dust Belt: 

Fact or Fiction. J. Geophys. Res. (to be published). 

2. Private discussions with personnel at Langley Research Center, Sept. 1965. 

3. Nilsson, C. : Some Doubts about the Earth's Dust Cloud. Science (to be 
published). 
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reasonably consistent with estimates based on penetration experiments and 
zodiacal light studies. However, there is still the unexplained fact that collection 
experiments flown on satellites as well as rocket probes have returned a very 
large number of particles, some of which were collected undamaged. 

It is clear that there are still many unanswered questions that must be 
resolved before all of the above measurements can be properly interpreted. 

Since there is some doubt concerning the measurements that imply the large 
concentration of non-penetrating particles, these data are not used in the present 
analysis. If these measurements can be verified, and the existence of a near- 
earth dust belt of extremely tenuous particles is demonstrated, the mass distri- 
bution obtained in this present analysis must be modified to include this com- 
ponent. Several challenging problems concerning the origin and nature of such 
particles must then be solved. It is not considered likely that such a finding 
will alter the penetration frequency predictions obtained by the present study 
for thicknesses greater than 0. 01 mm aluminum. 


CONCLUSIONS 


A penetrating flux model has been developed by essentially extrapolating 
the Pegasus data to larger and smaller thickness using slopes d log N/d log M 
determined from ground based and other satellite penetration measurements. 

The major source of uncertainty lies in the functional relations of mass to 
penetration and of mass to peak photographic magnitude. A secondary source 
lies in the freedom of choices in the representation of these slopes as a function 
of cumulative flux. This freedom can be restricted substantially by additional 
radar or photographic measurements of this log derivative from N = 10 4 to 10 e 
per m 2 day steradian, i. e. , between Elford's data and the Super-Schmidt data. 
Such measurements can certainly be made and perhaps already exist, but have 
not been reported in a form amenable to this analysis. As for the uncertainties 
in the mass relation, the problem arises because the empirical values for the 
mass exponent differ from those embodied in the theoretical models. Therefore, 
some uncertainty exists until the reasons for these differences are understood. 

A "best fit" cumulative mass flux has been found to within a constant of 
integration without having to specify the mass involved in any observation. A 
constant of integration was calculated for each measurement by estimating a 
value for the mass (having average velocity, average density, and normal in- 
cidence) characteristic to the observation. The average of these constants was 
used in the "best fit" mass flux. 
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The mass distribution function was obtained from the cumulative mass 
flux and is found to agree remarkably well with the interplanetary dust distri- 
bution derived by Van de Hulst from interpretation of photometric measurements 
of F-corona and zodiacal light. Also it was found that interplanetary measure- 
ments of dust by Alexander on Mariner IV were in order of magnitude agreement. 
These results suggest that the environment is not substantially different in 
interplanetary space. Also the fact that the slopes determined from satellite 
penetration data appear to rapidly approach zero as N increases, suggests an 
upper limit on the number of meteoroid encounters which appears to be on the 
order of 3/m 2 day regardless of size. This low value is consistent with the fact 
that only one discernable impact crater was found on a Gemini window. This 
also suggests that erosion from meteoroids is insignificant and that puncture 
rate in thin materials is far less than earlier estimates. This low value for the 
maximum encounter frequency demands that a meaningful meteoroid experiment 
must expose an area on the order of a m 2 to collect a good statistical sample in 
a reasonable time. 

There is a possibility that there exists, in addition to the family of 
meteoroids treated in this study, a very large population of particles concentrated 
near-earth that for some reason are not able to penetrate materials as thin as 
0. 01 mm aluminum. How such a concentration could exist and how the particles 
in this family could have such low penetrating power are not understood. Until 
the existence of this family is more conclusively demonstrated, it will not be con- 
sidered. Even if this family does exist, the low penetrating ability of its con- 
sistent particles must be such that it will not contribute significantly to the 
meteoroid hazard estimates presented in this paper. 
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